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Center. 
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Dr. J. J. Wortman of the Institute and Drs. L. K. Monteith and 
J. C. Mulligan of North Carolina State University made major contri- 
butions to the study. 

We also wish to acknowledge our appreciation for the technical 
advice and loan of some equipment from the Environmental Protection 


Agency. 



ABSTRACT 


A metal-oxide- silicon (MOS) capacitor-type particulate sensor was 
evaluated for use in atmospheric measurements. An accelerator system was 
designed and tested for the purpose of providing the necessary energy to 
trigger the MOS-type sensor. The accelerator system and the MOS sensor 
were characterized as a function of particle size and velocity. Diamond 
particles were used as particulate sources in laboratory tests. 

Preliminary tests were performed in which the detector was mounted 
on an aircraft and flown in the vicinity of coal-fired electric generating 
plants. 


iii 



STUDY OF THE USE OF METAL-OXIDE- SILICON (MOS) DEVICES 


FOR PARTICULATE DETECTION AND MONITORING 
IN THE EARTH'S ATMOSPHERE 

Table of Contents 

Preface ii 

Abstract iii 

List of Figures v 

List of Tables vii 

1.0 INTRODUCTION 1 

2.0 CALIBRATION PARTICLE SELECTION , 3 

2.1 PARTICLE DISPERSION METHOD 3 

2.2 PARTICLE SUSPENSION PREPARATION ... 7 

3.0 DIAMOND PARTICLE DISTRIBUTION SYSTEM ... 10 

4.0 CLIMET PARTICLE COUNTER 12 

5.0 MOS SENSOR FABRICATION 16 

6.0 VACUUM NOZZLE PARTICLE ACCELERATOR 20 

7.0 MOS ELECTRONICS . . 24 

8.0 DATA- GATHERING METHODS 28 

8.1 THRESHOLD DATA 30 

8.2 VARIABLE BIAS DATA 34 

8.3 VARIABLE PRESSURE DATA 34 

8.4 IN-FLIGHT SENSOR EVALUATION 38 

9.0 ANALYSIS 44 

10.0 DISCUSSION 50 

11.0 APPLICATION 56 

APPENDIX A: Gas and Particle Flow Characteristics 

for Suction Type Particle Accelerator 

used with MOS Detector 58 


List of Figures 


1 Size Distribution of 1-y Diamond Particles 4 

2 Size Distribution of 3-y Diamond Particles 5 

3 Size Distribution of 6-y Diamond Particles 6 

4 Schematic Representation of Particle Distribution System. ... 11 

5 Photograph of the Climet Particle Counter 14 

6 Block Diagram of the Climet Particle Counter 15 

7 Cross Section of MOS sensor 17 

8 Photomicrograph of MOS Sensor 19 

9 Modified Vacuum Nozzle Particle Accelerator 21 

10 Photograph of Vacuum Nozzle Particle Accelerator 22 

11 Sampling Characteristics of Vacuum Nozzle Particle 

Accelerator 23 

12 MOS Discharge Pulse 25 

13 Schematic Diagram of the MOS Readout Circuit 26 

14 A Schematic Representation of the Calibration System 29 

15 Particle Velocity as a Function of Particle Accelerator 

Nozzle Pressure (3 y, 9 y, and 15 y Diameter Particles) ... 31 

16 Particle Diameter-Mass for Diamond Microparticles 32 

17 Particle Mass-Velocity Threshold for Initiating MOS 

Capacitor Discharge 33 

18 Dependence of Velocity Threshold on Applied Voltage 35 

19 MOS Counting Efficiency Relative to Climet 6-Micron Particle 

Distribution (1200 A Si02) 36 

20 Photomicrograph of Impacted MOS Sensor 37 

21 Illustration of Geometrical Considerations to Determine 

Counting Efficiency 39 

22 Schematic of MOS Sensor Counting Circuit Used in Flight 

Evaluation 40 


v 




List of Figures (Continued) 


23 Block Diagram of MOS Sensor Counting System 40 

24 Photograph of MOS Sensor on Wing Spar of Aircraft 43 

25 Illustration of the Particle Penetration Geometry 46 

A-l Nozzle Flow and Velocity versus Nozzle Pressure 67 

A-2 Nozzle and Straight Tube Combination 69 

A-3 Particle Velocities as Functions of Nozzle Pressure 76' 


vi 



List of Tables 


1 Particle Counting Statistics and Accuracy of Climet 13 

2 Particle Penetration Threshold for Various Silicon 
Oxide Thicknesses with Electric Field of 2.2 x 10^ 

volts/cm 48 


vii 



1.0 INTRODUCTION 


The primary objective of this research effort was to construct, test, 
and calibrate a sensor to count microparticles from roughly 1 to 10 microns. 
Typical applications were considered to be the sensing of particulate 
matter from solid propellant rocket exhausts (alumina), particulate matter 
in the upper atmosphere of concern from supersonic transports, , and emissions 
from stationary sources such as power plants and from dusty environments 
such as rock quarries and coal mines . 

The microparticle sensor is a capacitor- type sensor where the 
impacting particle initiates a temporary short through the dielectric of 
the capacitor. The large current density through the shorted region of 
the capacitor raises the temperature in the region about the impact so 
that the material (dielectric, particle, electrode) in this region is 
vaporized. With the short circuit removed, the capacitor recharges. 

From penetration theory one expects that the microparticles at subsonic 
velocities will cause a discharge upon penetration to some fraction of 
their diameters. To insure adequate penetration to initiate an electrical 
discharge, the dielectric of the capacitor was chosen as thin as practical 
with modem silicon-silicon oxide technology. Silicon technology was 
chosen due to the sophistication and understanding available to achieve a 
very thin, reliable, and reproducible dielectric. Typically, the final 

O 

capacitor is approximately 1000 A, thermally grown silicon oxide film on 

O 

low resistivity (less than 1 ohm-cm) silicon wafer with a 500 A aluminum 
electrode which is impacted by the microparticle. 

The major problems encountered during this research effort were 
(1) obtaining a well-defined and reproducible particle source, (2) reducing 



unwanted and random sources of particles from the system, and (3) accelerat- 
ing particles to a velocity necessary to initiate a discharge. Much of the 
early effort was devoted to development of particle selection techniques. 
Well-defined polystyrene microparticles and organic particles were used; 
however, they would not initiate a reliable discharge. Our experience has 
shown that either density or hardness near that for alumina and diamond 
will initiate discharges. However, the major effort and probably the most 
significant result is the design of a microparticle dispensing system which 
provides a range of small particles with a degree of reproducibility 
unobtainable prior to this effort. Only through detailed design and 
analysis could the influence of factors such as unwanted microparticles 
throughout the system be identified. 

Finally, it should be pointed out that the sensor designed and tested 

in this research effort will count microparticles with a diameter greater 

3 

than 2 microns and a density near 3.5 gm/cm . The capacitor- type units 
were compared with a commercially available unit for counting microparticles, 
and the counts compared favorably and well within experimental error. 
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2.0 CALIBRATION PARTICLE SELECTION 


Diamond particles were chosen for the calibration source mainly for 
two reasons: (1) their hardness properties and (2) their narrow size 

distribution. Through previous experience with polystyrene, organic 
crystals, and alumina microparticles, the hardness of the particle is 
an important factor in obtaining reliable discharges of the impacted MOS 
capacitor. The narrow size distribution property was important in order 
to obtain a meaningful calibration of the MOS sensors for different size 
particles . 

The diamond particles were obtained from the Geoscience Instruments 
Corporation and were purchased in the following size ranges: 0.25 y, 

1.0 y, 3.0 y, 6.0 y, 9.0 y, and 15.0 y. Shown below is a size distribu- 
tion table for these particles. 


Average size 

0.25 y 

l y 

3 y 

6 y 

9 y 

15 y 

Distribution 

(0.1-5) 

(0.5-2) 

(2-4) 

(4-8) 

(8-12) 

(12-20) 


Also shown in Figures 1, 2, and 3 are distribution curves for the 1 y, 
3 y, and 6 y particles, respectively. 

2.1 PARTICLE DISPERSION METHOD 

The method chosen for dispersion of the diamond particles was to 
suspend the particles in isopropyl alcohol and to spray them into a 
sampling chamber with an atomizer. The suspension method was chosen 
because of the difficulty encountered in dispersing dry particles due 
to coagulation. Also, this method allowed much more accurate control 
of the number and uniformity of the dispersed particles. 


3 












Percent by weight 



Figure 2. Size Distribution of 3-y Diamond Particles 




Particle size in microns 


Figure 3. Size Distribution of 6-y Diamond Particles 





2.2 PARTICLE SUSPENSION PREPARATION 


Isopropyl alcohol was used for the suspension fluid because it is a 
volatile fluid with a comparatively high viscosity. As a first step in 
preparing the diamond particle suspension fluid, a particle background 
level in the isopropyl alcohol was determined by spraying the fluid into 
a sampling chamber with an atomizer driven by filtered nitrogen and 
sampling with a Climet Model 250 particle counter. 

Although the particle size distribution was very narrow, a sizeable 
number of smaller than nominal size particles in a given range were present. 
The method used for separating these smaller particles from the particles 
in the size range of interest was as follows. The settling velocity in 
isopropyl alcohol was calculated for the particles in the size range of 
interest. The particles were weighed out and placed in a glass container 
with a known volume of isopropyl alcohol. These were thoroughly mixed 
until they were uniformly distributed in the solution. The time for the 
particles of interest to settle to the bottom of the container was deter- 
mined using the height of the saturated fluid in the container and the 
calculated settling velocity of these particles. After the particles of 
interest had settled to the bottom of the container, the fluid containing 
the still suspended smaller particles was carefully siphoned into a 
graduated cylinder. This volume of fluid (isopropyl alcohol) was replaced 
by clean fluid. This procedure was repeated once in the case of the 6 p 
and 9 p particles and twice for the 3 p particles. This was necessary in 
order to obtain a solution of particles which would contain mainly the 
particles in the size range of interest. Shown below is the method used 
for determining the settling velocity and the settling time for the parti- 
cles of interest. Also shown is a table containing the calculated settling 
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velocities of various sized diamond particles in isopropyl alcohol. The 
settling velocity, V , is given by 


V 

s 



18n 


where p = density of particle 

g = acceleration due to gravity 
d = diameter of particle 
n = viscosity of suspension fluid. 


In air, 1-y spherical particles of unit density have a settling velocity 


_3 

V = 3.49 x 10 cm/sec . 
s 


In isopropyl alcohol, the settling velocity of the diamond particles may 
be obtained by introducing the ratio of the viscosity of isopropyl 
alcohol to air and the particle density where 


Viscosity of isopropyl 
Viscosity of air 


2 x 1Q~ 2 poise = 1>Q99 x 1q -2 
1.82 x 10 poise 


and 


Density of diamond = 3.5 , 


yielding 
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V 

s 


(3.5) (d^) (3.49 x 10 cm/sec)(6Q sec/min) 


1.099 x 10 


-2 


= (d^)6.67 x 10 ^ cm/min . 


Particle size (y) 

0.5 y 

1 M 

3 y 

6 y 

9 y 

15 y 

V (cm/min) 
s 

0.0016 

0.0066 

0.06 

0.24 

0.54 

1.5 


„ . . . height of fluid (cm) 

Settling time = — — s - 

s 
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3.0 DIAMOND PARTICLE DISTRIBUTION SYSTEM 


Due to the more rapid settling velocity and aerodynamic properties of 
the diamond particles, the Environmental Research Corporation fluid atomi- 
zation generator, originally planned to be used to disperse the particles, 
would not efficiently disperse the large particles. Because of this, a 
particle dispersing system was designed which overcame these problems 
and would dispense the large diamond particles. An illustration of this 
system is shown in Figure 4. The isopropyl alcohol solution containing the 
suspended diamond particles is stirred by a magnetic stirrer to maintain a 
uniform suspension of particles and is then fed by a combination of gravity 
and venturi action to the nozzle in the sampling chamber. By using the 
gravity feed technique, one is assured of getting the large particles into 
the chamber. The venturi-type atomization nozzle was fed by nitrogen fil- 
tered by a 0.45-y filter at 30 psi. The nitrogen flow rate through the 
nozzle was measured by a wet test meter to be 21 liters/minute. The high 
velocity stream of nitrogen through the venturi nozzle caused the isopropyl 
alcohol solution containing the suspended diamond particles to be sprayed 
into the chamber in a fine mist. The alcohol evaporated before striking 
the opposite wall of the sampling chamber, allowing the particles to settle 
by gravitational action to the bottom of the sampling chamber where they 
were simultaneously sampled by the Climet Model 250 particle counter and 
the MOS vacuum nozzle particle accelerator. An additional 0.45-y filtered 
nitrogen source of 18 liters/minute was added to the chamber so that the 
total nitrogen input was approximately a factor of two greater than the 
combined sampling rates of the Climet and the vacuum nozzle particle 
accelerator. This technique insures that a slight positive pressure re- 
mains in the sampling chamber, thus eliminating any effect from the 
surrounding atmosphere. 
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Figure 4. Schematic Representation of Particle Distribution System 



4.0 CLIMET PARTICLE COUNTER 


The Climet Model CI-250 particle counter was used as a reference 

during the calibration of the MOS particulate sensors. The Climet particle 

counter is a portable, solid-state, optical instrument operating on the 

light-scattering principle. The air or gas containing the particles to be 

sampled is pulled into the instrument by a vacuum pump at a rate of 
3 

14 ft /hr. The large sampling rate is to insure that all the particles 

in the sample will enter the instrument. The sampling air stream is 

3 

divided so that 1 ft /hr of the sample is pulled through the optical 
counting chamber. Any particle which enters the counting chamber causes 
the light to be scattered. The scattered light pulse, which has an 
amplitude proportional to the size of the particle, passes through a lens 
system and strikes a silicon photocell. The electrical pulse from the 
photocell is fed to a comparator circuit which is controlled by the 
particle size range switch on the instrument. Any particles which are 
equal to or larger than the particle size range switch setting produce 
pulses which trigger the comparator circuit which in turn drives the 
digital counter. The counting time is controlled by selecting a 36-second, 
a 360-second, or a manual mode of operation. The accuracy of the instru- 
ment is illustrated in Table 1. A picture of the instrument and a block 
diagram are shown in Figures 5 and 6, respectively. 

The instrument was modified for use in this study. It was found that 
the large particles (75 p) would not reach the counting chamber through 
the normal plumbing. A plastic pipe was inserted directly into the chamber 
and positioned in a vertical manner so that particles would fall directly 
from the source to the chamber. 
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Table 1. Particle Counting Statistics and Accuracy of Climet 


— ; 1 

Number of 
particles 
counted 

Percent error 

Standard 
deviation 
of counts, 3 

50% proba- 
bility 

95% proba- 
bility 

99.7 % proba- 
bility 

99,999 

±0.213 

• ±0.63 

±0.95 

316 

10,000 

±0.67 

±2.0 

±3.0 

100 

1,000 

±2.13 

±6.3 

+9.5 

31.6 

100 

±6.7 

±20 

±30 

10 

10 

±21.5 

+65 
. -60 

+ 100 
- 75 

3.16 

1 

+ 200 
- 70 

+ 300 
| - 90 

+500 

-100 

1 
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Figure 6. Block Diagram of the Climet Particle Counter. 




5.0 MOS SENSOR FABRICATION 


The MOS particle sensors are basically MOS capacitors fabricated by 
using standard integrated circuit processes. The basic structure of the 
MOS sensor is a silicon substrate acting as one "capacitor plate," a 
thermally grown Si 02 layer acting as a dielectric, and an evaporated 
aluminum film acting as the second capacitor plate. A cross-section of 
the unit is shown in Figure 7. The silicon substrate material is p-type 
0.1 p-cm, 1-1-1 oriented, 1-1/4- inch diameter wafers polished on one 
side. Following is a description of the MOS sensor fabrication process. 

The first step in the fabrication process was to prepare the silicon 
substrate and grow the Si02 film. The silicon substrates were cleaned by 
employing standard silicon cleaning procedures which included heating the 
substrate in solvents and acids. The oxidation step was carried out by 
loading the substrates into a quartz oxidation tube of an oxidation fur- 
nace which was at 1100°C, with 2000 cc of dry 0^ flowing through it. The 

O 

time of oxidation was varied in order to produce Si02 films of 700 A, 

o o o o 

1000 A, 1200 A, 1500 A, and 1800 A. 

The next step was to anneal the substrates containing the Si02 layers 
at 200°C in a vacuum of 25 in. Hg for 30 minutes. The wafers were then 
placed in the aluminum evaporation system on evaporation masks which con- 
tained a 50 x 50 array of 15-mil diameter holes on 20-mil centers. A 

O 

standard aluminum evaporation cycle was carried out which produced 500 A 
thick, 15-mil diameter dots of aluminum on the Si02 surface. 

The next step in the fabrication sequence was to etch the Si02 from 
the lapped back side of the wafer by positioning the wafers, aluminum 
side down, and placing crystals of ammonium bifluoride (NH^F’HF) on the 
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- 500 A A1 (front electrode) 



SiO^ (dielectric) 


Si (rear electrode) 


Figure 7. Cross Section of MOS Sensor. 




lapped wafer surface. A few drops of D.I. 1^0 were placed on the crystals, 
and when the SiC^ was etched away, they were rinsed in D.I. H^O. 

O 

The next step was to evaporate 500 A of aluminum on the etched back 
surface of the silicon wafers. After this step, the silicon wafers were 
scribed into 100 x 100 mil square chips containing approximately a 5 x 5 
array of 15-mil diameter dots. These 100-mil chips were eutectically 
bonded to 12 pin gold-plated TO-5 headers by scrubbing the chip in a 
random motion on the T0-5 header at 450°C. The eutectic bond formed was 
mechanically strong and produced a good electrical contact between the 
silicon chip and the TO-5 header. 

The next step was to bond 1-mil gold leads from eleven of the 15-mil 
aluminum dots to the posts on the TO-5 header. The twelfth post was 
bonded to the base of the T0-5 header and served as a common contact for 
the eleven 15-mil diameter MOS capacitors. A photomicrograph of the com- 
pleted MOS particle sensor is shown in Figure 8. 
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Figure 8. Photomicrograph of MOS Sensor. 
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6.0 VACUUM NOZZLE PARTICLE ACCELERATOR 


The purpose of the vacuum nozzle particle accelerator is to accel- 
erate the particles to a velocity sufficient to produce counting on the 
MOS particle sensor. A sketch and a photograph of the accelerator are 
shown in Figures 9 and 10, respectively. The vacuum nozzle particle 
accelerator was chosen, as a means of obtaining the necessary particle 
acceleration, because of its simplicity of design and compact size. 

The design criteria are given in the Appendix. The particle acceleration 
tube was added to the vacuum nozzle particle accelerator in order to pro- 
vide sufficient distance, and thus time, for the 6 p and larger diamond 
particles to reach a velocity sufficient to produce counting on the MOS 
sensor. The T0-5 MOS particle sensor was mounted into a standard 12-pin 
TO-5 I-C socket which positioned and aligned the sensor in front of the 
particle acceleration nozzle. The alignment of the MOS particle sensor 
was found to be a very important factor during the design of the vacuum 
nozzle particle sensor. The 0.5- inch nozzle-to-MOS sensor distance was 
chosen in order to prevent small, soft particles of insufficient mass from 
producing counts and from accumulating on the MOS sensor surface, thus 
lowering the counting efficiency of the sensor. In operation, a vacuum 
pump was connected to the vacuum port on the vacuum nozzle particle 
accelerator, and the vacuum nozzle pressure was measured at the vacuum 
port. A plot of the sampling rate versus the vacuum nozzle pressure is 
shown in Figure 11. 


20 




o <u 

•H A! 


4-1 CO 


U 4-1 


CO C 


Pw -H 


21 


Figure 9. Modified Vacuum Nozzle Particle Accelerator. 






Figure 10. Photograph of Vacuum Nozzle Particle Accelerator. 
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Figure 




7.0 MOS ELECTRONICS 


The MOS particle sensor is, in effect, a MOS capacitor having the 
silicon substrate as one plate, the Si 02 film as the dielectric, and a 

O 

500 A film of evaporated aluminum as the other plate. A DC bias, sup- 
plied by a battery, is placed across the MOS capacitor through a 1000 ft 
resistor. When a particle strikes the MOS sensor with sufficient momentum 
to cause damage of sufficient magnitude in the dielectric, the capacitor 
discharges through the dielectric. The capacitors are connected in 
parallel to obtain large counting areas with the added benefit that the 
parallel connection provides added charge, hence current density, to burn 
away the aluminum in the damaged area. The time for this clearing opera- 
tion is on the order of nanoseconds as shown by the rise time in Figure 12. 
The battery then recharges the capacitor through the 1000 ft resistor and 
thus reestablishes the electric field across the capacitor as shown by the 
fall time of the pulse in Figure 12 . A portion of the recharge current is 
fed to a Fairchild pA-710 comparator circuit. Since there is a variation 
in the amplitude of the recharge pulses, which depends upon the applied 
bias and to some extent upon the impact event, this circuit is necessary 
in order to supply constant amplitude pulses to the HP 25456 electron 
counter. For most instances, the amplitude of the discharge pulse is 
approximately 25 percent of the applied bias. A schematic diagram of the 
MOS readout circuit is shown in Figure 13. Shown below is a sample cal- 
culation of the capacitance and maximum counting rate of the MOS sensors. 

C(pf) = 8.85 x 10 -2 (!) 

d 
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where K = 3.2 for SiO^ 

d = dielectric thickness in centimeters 

. . 2 
A = area m cm . 

-3 2 

For a 15-mil diameter MOS unit, the area is 1.14 x 10 cm . With 
° -5 

a thickness of 1000 A (10 cm), the 15-mil diameter MOS capacitor has a 
capacitance of 32.3 picofarads. Assuming that eleven units are connected 
in parallel, as in most instances, the total capacitance, C^, is 355 pico- 
farads. Thus the RC time constant for recharging the impacted sensor is 

RC = (1 x 10 3 «) (355.1 x 10 -12 f) = 0.355 x 10~ 6 sec . 

Allowing one RC time constant for each counting event yields 

11 6 
f = — = —7 = 2.8 x 10 maximum counts/sec . 

P 0.355 x 10 b sec 

Below is a table showing the effect of SiO^ thickness on the MOS 
capacitance and maximum counting rate. 


Si 02 thickness 

C 

C T 

RC 

f 

700 A 

46.1 pf 

507.1 pf 

0.51 ysec 

1.96 x 10 6 

O 

1000 A 

32.3 pf 

355.3 pf 

0.36 ysec 

2.77 x 10 6 

1200 A 

26.9 pf 

295.9 pf 

0.29 ysec 

3.39 x 10 6 
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8.0 DATA-GATHERING METHODS 


The calibration procedure was carried out by simultaneously sampling 
from the particle generation system with the vacuum nozzle particle accel- 
erator and the Climet Model CI-250 particle counter acting as a reference. 
The calibration setup is illustrated in Figure 14. 

At the beginning of each set of runs, the background of the system 
was checked by using isopropyl alcohol in the particle dispensing system 
and sampling simultaneously with the Climet particle counter and the MOS 
sensor in the vacuum nozzle particle accelerator. In general, the Climet 
would indicate a background count which the MOS would not. The back- 
ground for the Climet could then be subtracted from the Climet counts 
obtained from the runs made with a particle-saturated solution to obtain 
the correct number of counts to be used as a standard for comparison with 
the MOS counts. Each time a run with different size particles was made, 
the dispensing system was completely cleaned and the background was 
checked . 

The three basic types of data taken were (1) threshold data, 

(2) variable pressure data, and (3) variable bias data. In addition, a 
preliminary evaluation as an environmental monitoring instrument was con- 
ducted by mounting a MOS unit aboard an aircraft and flying through the 
emissions from power plant smoke stacks in Virginia and North Carolina. 

The variables examined during these series of tests were dielectric (SiO^) 
thickness, particle size, particle velocity, sampling rate, and MOS bias 

O O 

voltage. The three dielectric thicknesses investigated were 700 A, 1000 A, 

O 

and 1200 A; and the diamond particle sizes investigated were 1 p, 3 p, 6 p, 
and 9 p. 
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Figure 14. A Schematic Representation of the Calibration System. 




8.1 THRESHOLD DATA 


The threshold data were obtained to determine a mass-velocity rela- 
tionship for particle counting by a MOS sensor with a given dielectric 

O 

thickness. Threshold data were obtained for the 700, 1000, and 1200 A 
silicon oxide thicknesses and for 3-, 6-, and 9-micron diamond particles. 
With the MOS sensor biased at a given voltage, the threshold for counting 
was determined by starting with the MOS vacuum nozzle particle accelerator 
velocity at zero and then slowly increasing the nozzle pressure, thus the 
particle velocity, until the MOS sensor began to count. For each set of 
conditions — bias, particle size, and Si 02 thickness — the counting thresh- 
old nozzle pressure was obtained ten times and then averaged. The 
nozzle pressure was converted to a particle velocity using the theory in 
Appendix A. The graphical representation of the theory used is displayed 
in Figure 15 where the ratio of nozzle pressure to atmospheric pressure 
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Figure 15. Particle Velocity as a Function of Particle Accelerator Nozzle Pressure 
(3 y, 9 y, and 15 y Diameter Particles). 
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Figure 17. Particle Mass-Velocity Threshold for Initiating MOS Capacitor Discharge. 



8.2 VARIABLE BIAS DATA 


Variable bias data were taken to determine the effect of an electric 
field upon the threshold of counting for the MOS sensor. Observing the 
threshold for several bias voltages with each particle distribution and 
with each thickness of SiO^ resulted in the data included in Figure 18. 

The straight lines represent a theoretical interpretation of the data to 
be discussed in a later section. 

8.3 VARIABLE PRESSURE DATA 

The variable nozzle pressure data were taken to determine the count- 
ing efficiency of the MOS sensor relative to the Climet instrument. In 
addition, the particle size selectivity of the MOS can be inferred assum- 
ing a correlation between the instruments exists and assuming a consistent 
model for the operation of the MOS can be discerned. A typical set of 
data is shown in Figure 19. These data were obtained by varying the 
nozzle pressure and noting the MOS particle count. At least five separate 
observations at each nozzle pressure were averaged to obtain each data 
point shown in Figure 19. Each observation was for a 36-second sampling 
interval. For each sampling interval a Climet reading in either the 5- or 
2-micron range was obtained. For the particular range the instrument 
counts all particles about that setting. Since the volume of air sampled 
by the Climet and the MOS in the 36-second interval differed, a correc- 
tion factor to obtain a meaningful comparison was introduced. Basically, 
the average of the Climet counts was multiplied by the ratio of the volume 
of air sampled by the MOS to the volume of air sampled by the Climet . In 
addition, the MOS counts were corrected for the fact that the active area 
of the device is only a fraction of the area intercepted by the sampled 
air stream. An example of the area effect can be seen in Figure 20, where 
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Figure 20. Photomicrograph of Impacted MOS Sensor. 
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the number of 15-mil dots that show a preponderance of counts is limited 
to three or four units. Using simple geometric considerations as shown 
in Figure 21, the average of the MOS counts at each nozzle pressure was 
corrected. Through actual examinations of the surface of an impacted 
MOS sensor, it was determined that somewhere between three and seven of 
the 15-mil MOS units did most of the counting. Thus the MOS counts should 
be multiplied by a factor between approximately 2 to 4. Numerous observa- 
tions led to the use of the larger correction of 4. Using the above 
procedures certainly limits the accuracy of the results; however, quali- 
tative interpretation predicated on these results will be given in a later 
section. 

8.4 IN-FLIGHT SENSOR EVALUATION 

The purpose for conducting the in-flight tests of the MOS sensors 
was to determine the ability of these sensors to function as useful 
instruments for environmental monitoring. 

A Cessna 182 was chartered for use on these test flights. Two 
vacuum nozzle particle accelerators were mounted on the wing spar. The 
particle intake ports were aimed toward the rear of the plane to prevent 
dust from the runway from damaging the MOS sensor during takeoff. A 
triggering circuit, illustrated in Figure 22, was designed and used 
with the normal MOS sensor circuitry in order to trigger a strip chart 
recorder which acted as an event counter. The AC power for the vacuum 
pump and the electronic equipment was supplied by two 12-volt wet cell 
batteries and a 500-watt inverter. The MOS sensors used for these tests 

O 

were the 4><4 array, 15-mil diameter dot type with 500 A aluminum con- 

O 

tacts. The Si02 dielectric thicknesses were 700 and 1200 A. A block 
diagram of the equipment used for these tests is shown in Figure 23. 
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vacuum nozzle Dimensions in mils 


Figure 



Illustration of Geometrical Considerations 
to Determine Counting Efficiency. 
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Figure 22. Schematic of MOS Sensor Counting Circuit Used 
in Flight Evaluation. 



Figure 23. Block Diagram of MOS Sensor Counting System. 
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Since the MOS particle sensors had individual lead-in cables and vacuum 

O O 

lines, either the 700 A or the 1200 A sensor could be easily switched 
into the counting circuitry. 

The MOS particle sensors were in operation during the flight from 
the Raleigh- Durham Airport to Langley Field and on the return trip to 
Raleigh-Durham. Three flights were made from the Langley Field over the 
Yorktown area. The flights were made at altitudes ranging from 500 to 
1100 feet. The particle counting rate was largest at 500 feet and 
dropped off with increasing altitude. Particle counts were recorded 
more than 2 miles downwind from a coal-fired power plant. The visible 
plume from the source was approximately 1000 feet high; however, the 
larger counting rate at 500 feet indicated that the larger particles had 

O 

settled out of the plume. Fourteen counts were recorded with a 1200 A 
MOS sensor as the stack was approached from downwind at 1000 feet alti- 
tude, with the highest fate of counting being approximately 0.5 mile from 
the stack. At an altitude of 800 feet, an identical approach was made 

O 

toward the stack, and a total of 61 counts were recorded. When the 1200 A 
MOS sensor was flown at 500 feet at right angles to the plume, 523 counts 

O 

were recorded. The counting rate at 500 feet with a 700 A MOS sensor was 
of sufficient magnitude to destroy the sensor after a few seconds of 
counting. 

A light rain occurred between the second and third flights, and the 
particle counting rate dropped off by approximately a factor of ten. The 
rain droplets apparently condensed on the particles and thus cleaned the 
air as the rain fell. No evidence of corrosion of the metal (aluminum) 
contacts was observed, indicating that this anticipated problem did not 
materialize. 
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A flight over a North Carolina coal-fired electric power plant 

O O 

yielded data with both the 700 A and 1200 A MOS sensors that were very 
similar in magnitude and altitude conditions to those recorded in 
Virginia. A photograph showing the mounting arrangement on the wing spar 
of the aircraft is shown in Figure 24. 
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Figure 24. Photograph of MOS 
Sensor on Wing Spar of Air- 


9.0 ANALYSIS 


To facilitate calibration of other MOS units, a consistent model which 
relates the experimental observations in the present effort will be 
presented. The threshold data in Figure 17 suggest a rather classical 
penetration relationship where 


X = k m^ v 11 (2) 

P 

where is the thickness penetrated by a particle of mass m and velocity v. 
A reasonable assumption to explore is that penetration through the aluminum 
electrode and the dielectric is required to initiate an electrical dis- 
charge. However, two observations are not included in such a simple model. 
First, it has been observed that the threshold is dependent upon the 
applied bias. As voltage is increased, the threshold velocity for a given 
particle mass decreases. For typical data in Figure 17 the particle energy 
is near 10 ^ joules while the energy stored on the capacitor is near 
10 ^ joules. Thus the electrical energy is between 100 and 1000 times the 
kinetic energy of the particle. Thus the added particle energy as the 
velocity is increased does not replace the reduction in electrical energy 
with reduced applied voltage. Secondly, a rather unique feature of the 
operation of the MOS units is that for applied voltages less than 7 volts, 
a particle impact results in a short which will not clear. This lower 
limit of operation is independent of oxide thickness at least for the units 
used in this study. 

Another factor which will be useful in characterizing the operation 
of the MOS sensor is the inherent electrical field strength for the silicon 
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oxide. In numerous observations the units would .spontaneously discharge 

when the electrical field was between 3 and 4 x 10° volts/cm. Detailed 

6 7 

studies of electrical breakdown in SiC^ report values from 10° to 10 / 
volts/cm. At the lower values it is generally recognized that small 
defects in oxide are responsible for the discharges. The electrical 
discharges are indistinguishable from the microparticle- initiated dis- 
charges and clearly must be avoided. Thus the range of operation is 

£ 

bounded on the lower end by the non-clearing events at about 10° volts /cm 
and on the upper end by the spontaneous electrical discharges at about 
3 x 10^ volts/cm. 

Although the present understanding evolved in a rather circuitous 
manner, the penetration- type relationship involving mass-velocity sug- 
gested that a well defined penetration was required. The voltage dependence 
of threshold suggested that penetration to a depth where the electrical 
field in the unpenetrated oxide exceeded the field strength for breakdown 
should be considered. A simple illustration of the geometry of the pene- 
tration is shown in Figure 25. As the particle penetrates to depth X^, 

' f \ 

it is assumed that aluminum is carried along to provide electrical cont-i-“. 
nuity such that the applied voltage appears across X^. When the depth of 
penetration is such that the field in X^ exceeds the field strength of the 
oxide, a discharge will occur. A rather tenuous assumption of a field 
strength near 10^ volts/cm could be used to proceed. However, the data in 
Figure 18 can be used to infer a value for the field strength. Assuming a 
penetration relationship for threshold of particle- initiated discharges 
where 

X = k m^ v 11 (3) 

P 


A 
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Figure 25. Illustration of the Particle Penetration 
Geometry. 
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( 4 ) 


and a field strength for discharge events (E ) where 


E = — — 
d X 3 ’ 


the dependence of threshold velocity upon applied voltage is given by 


V = - E, k m £ v n + E , (X- + X.) . 
d d 1 Z 


(5) 


Without attempting a detailed fit to the data, a value of = 4 x 10 
volts/cm was used. Partial justification for the choice was the observa- 
tions of spontaneous discharge near 4 x 10^ volts/cm. While field 
strengths near 10^ volts/cm have been reported, the impact event may 
damage the unpenetrated region resulting in a lower field strength. In 
any event, the value of 4 x 10 ^ volts/cm provides a reasonable lower limit. 
After assuming a value for E^, it is straightforward to find X^ for the 
mass-velocity relationship in Figure 17. Since the data were taken at 
equal applied electric fields of 2.2 x 10 volts/cm, the ratio X^/X^ is 
given by 


*2 = 4 x IQ 6 

X 3 2.2 x 10 6 


( 6 ) 


= 1.8 . 


Therefore, 


X = X, + X. - X- 
p 1 2 3 


- x i + x 2 - rh = x i + °- 4 < x 2 


(7) 
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The thickness of the aluminum electrode was 500 A in every case. The pene- 
tration threshold with an applied electric field of 2.2 x 10^ volts/cm for 
the three Si02 thicknesses is given in Table 2. 


Table 2. Particle Penetration Threshold for 
Various Silicon Oxide Thicknesses 
with Electric Field of 2.2 x 10^ 
volts/cm 


SiO„ thickness 
X 2 

Penetration thickness 
X 

P 

O 


700 A 

810 

1000 A 

940 

1200 A 

1030 


Using the penetration thicknesses in Table 2 to obtain a fit to the data 
in Figure 17 results in 


X = 8.5 x 10 4 m°’ 54 v (8) 

P 

where is in centimeters, m is in grams, and v is in centimeters /second. 
This linear dependence of penetration on velocity has been previously 
observed for impacts at low velocities. ^ 

Using the results in Equation 8 and a value of E^ = 4 x 10 volts/cm, 
a set of theoretical curves have been drawn in Figure 18. Some adjustment 
in E^ may result in a better fit; however, the adjustment would appear to 
be less than 10 percent based on the data for 6-micron particles.’ Howevhr, 
the 9-micron particle data are a bit high and suggest a slight increase. ' * 
It should be pointed out that the data were taken on a .number of units, 
and any variations in the field strength, from unit to unit will influence 
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the results. Therefore, some of the scatter in the data may result from 
variations in electrical properties. As mentioned earlier, the particle 
distributions are not sufficiently characterized to provide precise analysis 
therefore, the theoretical curves are considered reasonable under the 
circumstances . 
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10.0 DISCUSSION 


The calibration data in Section 8 and the analysis in Section 9 can 
be combined to provide an interpretation of the counting efficiency of the 
MOS relative to the Climet and to provide an estimate of the minimum 
particle diameter counted by the MOS. With a given bias on the MOS, the 
impacting particle must penetrate well into the oxide to initiate break- 
down. For the MOS particle accelerator analyzed in Appendix A, the 
velocity of the impacting particles is dependent upon the particle diameter. 
Thus at a given nozzle pressure, only particles of a certain diameter or 
greater will achieve the penetration threshold for discharging the MOS 
capacitor. The accelerator is essentially a mass selector in terms of 
threshold since an increase by a factor of ten in mass results in less than 
a factor of three change in velocity for a given nozzle pressure in 
Figure 15. If all particles impacted with identical velocities, the mass 
selectivity would be very precise. However, the threshold of penetration 
is approximately equal to the square root of the particle kinetic energy; 
thus a range of particles will have the same threshold of penetration for 
the MOS unit. Since the particle mass varies as the cube of the diameter, 
the threshold of penetration in terms of particle diameter will have a 
much narrower range. For a given nozzle pressure in the range of interest, 

0. 7 < (P„/P ) < 0.95, the particle velocity decreases approximately as the 

_> atm ~ ' 

inverse square root of the particle diameter, or 


V 

P 


'V 


d l/2 • 
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Thus the threshold penetration varies essentially as 


2 1/2 

X ^ (m v ) v d 
P P 


( 18 ) 


or directly proportional to particle diameter. This simply means that as 
the nozzle pressure is increased, the threshold will be determined by the 
largest diameter particles present in the sampled air volume. This 
assumes that the pressure is increased slowly enough that the relative 
occurrence of the particle diameters does not influence the measurement. 
Another way of stating this assumption is that over the observation time, 
all particles of various diameters in the sampled volume of air should 
impact a MOS capacitor. 

The counting efficiency observed for the 6-micron particle distribu- 
tion shown in Figure 19 can be explained in terms of the threshold of 
penetration dependence on particle diameter. As the nozzle pressure is 
increased, a threshold at 5 inches of mercury is due to the largest 
diameter particle in the particles sampled. By direct observation of par- 
ticles collected on a glass slide, this particle diameter is near 6 microns. 
The Climet reading for particle diameters greater than 5 microns is used as 
a standard, and the MOS counts 65 percent of the greater than 5 micron par- 
ticles at the threshold. As the nozzle pressure is increased, the MOS 
counts not only particles greater than 5 microns but also particles less 
than 5 microns in diameter, hence an efficiency greater than 100 percent 
since the 5-micron range on the Climet is used as the standard. Changing 
to the greater than 2 micron range on the Climet, an increase in nozzle 
pressure continues to increase the counting efficiency. This increase is 
due to the fact that the particles are distributed in diameter and as the 
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nozzle pressure is increased, smaller diameter particles reach the threshold 
of penetration and initiate a capacitor discharge. The fact that the MOS 
counted only 70 percent of the Climet counts may be due to either one of 
two possibilities: The precision of the data must surely be in question 

considering the uncertainties in the conversion factors used to provide a 
meaningful comparison between the Climet with the MOS, or there could be 
30 percent of the particles greater than 2 microns which do not achieve 
the threshold of penetration. Assuming the model is consistent, the 
smaller particles near 2 microns would constitute the 30 percent not 
counted. Using the data in Figure 15 and interpolating the velocity for 
a 2-micron particle at 19 inches of mercury indicates that the threshold 
of penetration is achieved; thus the 30-percent difference may simply be 
due to the lack of precision in correcting the data. 

The lower limit for counting with the MOS unit designed and employed 
for these calibration studies can be inferred. The airstream velocity 
is approximately 2900 cm/sec at 19 inches of mercury which is the maxi- 
mum achievable velocity for the present design. Assuming the particles 

O 

could be accelerated to this maximum velocity, the 700 A Si02 unit would 

O 

count all particles greater than 1.4 microns in diameter, the 1000 A Si02 
would count all particles greater than 1.5 microns in diameter, and the 

O 

1200 A Si02 would count all particles greater than 1.6 microns in diameter. 
Since the particles do not achieve airstream velocity, the lower limit is 
somewhat larger than the above values; however, smaller than 2 microns. 

To provide a test of these conclusions, a distribution of small particles 
which contains no particles greater than 2 microns was impacted against the 

O 

1000 A MOS units. The nozzle pressure was increased to 19 inches of mercury 
however, the MOS unit did not count. Climet readings were used to determine 


52 



the particle size distribution. To insure that the MOS counter would 
count, a distribution containing greater than 2 micron diameter was 
impacted. The MOS unit did respond to the distribution containing greater 
than 2-micron diameter particles. 

Regarding the observation that polystyrene particles do not discharge 
the MOS sensor, classical penetration theory would include the effect of 
the density of the particle and the density of the target material. In 
fact, two particles of density and impacting identical targets at 
identical velocities would result in ratios of penetration given by^^ 



Thus less dense particles penetrate less than the denser particles. Part 
of our observations may be simply explained by the density of the particles. 
However, in some instances with large organic crystals ( v 10 microns) which 
did not initiate discharges, the density of the particle does not account 
for the failure to initiate a discharge. These observations suggest that 
hardness must also be considered. We propose that brittle particles 
absorb some of the energy of the impact by shattering at impact, thereby 
reducing the energy available for the penetration. However, to validate 
the observation and the explanation will require further study. 

Finally, the shorting of the capacitors observed in this study when 
the applied bias is reduced to approximately 7 volts and threshold of 
velocity should be considered. For the different particles, they must 
achieve the penetration threshold regardless of particle size. At 
7 volts the threshold penetration is approaching complete penetration; 
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hence the discharge is not initiated by dielectric breakdown but by leak- 
age through a damaged dielectric. It is possible that the lack of a 
sharp rise time for the discharge precludes the rapid rise in local 
temperature, hence the lack of vaporizing the aluminum carried through 
by the particle. In any event, the unpenetrated dielectric which is 
assumed to be damaged and leaky rather than insulating increases by 30 

O O 

percent from a 700 A Si02 film to a 1200 A Si02 film. Considering the 
lack of precise control of particle velocity and size which impacts the 
unit with near threshold penetration, it is unlikely that a 30-percent 
variation when shorting is the ultimate observation can be established. 

In conclusion, it is unlikely that data can be obtained with sufficient 
precision to see the 30-percent variation. However, the range from 7 
to 8 volts as the applied bias which will result in shorted MOS units has 
been repeatedly observed. 
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11.0 APPLICATION 


In this report we have reported on the design, testing, and calibra- 
tion of a MOS particulate detector. Numerous variations in detector 
housing, in the MOS capacitor geometry, in the particle sampling tech- 
nique, and in the calibration procedures have resulted in the final design. 
We do not claim that the design has been optimized for every application; 
however, we consider the unit on which the report is based to be near an 
optimum design for detecting particles with a diameter between 2 and 10 
microns . 

Since the design of the particulate detector is to a large extent 
empirical, one can only conjecture about the design of a detector for 
particles greater than 10 microns and less than 2 microns. However, we 
are not aware of any fundamental reasons why a detector cannot be designed 
to operate outside the limits of 2 and 10 microns. The point is that for 
a detector to count particles with diameters from 2 to 10 microns, one 
should start with the unit described in this report. If the application 
is outside the 2 and 10 micron limits, from this report one can infer only 
a few factors pertinent to the design. 

The major factors to consider when applying the results included in 
this report are: 

(1) Geometric factors 

(a) The housing design should replicate the illustration in 
Figure 19. In particular, the dimensions of the particle 
accelerator tube and the nozzle to MOS sensor spacing 
should be maintained. 

(b) The MOS fabrication should follow the procedures in 
Section 5.0. Specifically, the aluminum thickness should 
be 500 A and the Si02 thickness should be 700, 1000, or 
1200 A. Aluminum dot area and spacing are critical; 
however, some variation can be considered since calibra- 
tion for efficiency will account for these variations. 
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(2) Particle Sampling 

The volume of air sampled will depend upon the nozzle pressure 
as shown in Figure 11. The nozzle pressure must be adjusted 
to insure that the particle you wish to detect achieves suffi- 
cient velocity to initiate a discharge. For example, if a 
9-micron particle is of interest, one finds in Figure 16 that 
the diamond microparticle has a mass of approximately 10“9 
grams. From Figure 17, the threshold velocity for a 10 - 9 gram 
particle impacting a 700 A oxide and 500 A aluminum electrode 
- is approximately 2,000 cm/sec. (This assumes a field of 
2.2 x 106 volts/cm in the oxide.) In Figure 15, a particle 
velocity of 2,000 cm/sec can be achieved with a nozzle pressure 
of 0.95 atmospheres. Thus the vacuum line should be reduced 
below this level for reliable counting. 

(3) Calibration 

(a) It has been inferred that particle hardness is important 
in the discharge process. To date most data have been 
obtained with alumina or diamond dust. For other dust 
environments, one should perform controlled laboratory 
tests with known (qualitatively) environments including 
the particles of interest. 

(b) The MOS capacitor may short after extended use. A short- 
ing test such as a current or voltage measurement across 
the MOS unit should be included in the counting circuit. 

A threshold of acceptable performance can be established 
on the basis of d.c. or average values which would pre- 
clude the observation of the transient discharge. 

(c) The counting rate should be less than 10^ counts per 
second. Even at this rate the MOS capacitors may degrade 
in a relatively short period of time. In most of this 
work we have considered 10 to 100 counts per second as 
acceptable counting rates which provide minutes to hours 
of continuous sampling without seriously degrading the 
response of the detector. 

(d) The threshold dependence on the voltage of the MOS capaci- 
tor will influence counting. One should be sure to use an 
applied voltage such that the field in the oxide is in 
excess of 2.2 x 10^ volts/cm. 

(e) Counting efficiency will depend upon the ratio of the 
active (electrically connected) MOS capacitor area impacted 
to the total MOS area impacted. Although we find the fac- 
tor of 4 (multiply counts by 4 to obtain actual count) to 
be a reasonable number for 15-mil dots on 20-mil centers, 
one should impact the sensor and then examine the impacted 
areas under a microscope to obtain assurance. Alignment 
errors can result in larger variations. One should remember 
that the entire area of the MOS array is the same order as 
the area of the nozzle outlet. Alignment must be carefully 
examined. 
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APPENDIX 


GAS AND PARTICLE FLOW CHARACTERISTICS 
FOR SUCTION TYPE PARTICLE ACCELERATOR 
USED WITH MOS DETECTOR 


By J. C. Mulligan 
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NOMENCLATURE 


A 

a 

b 

C. 


D 

d 

f 


w 


k 

K 

K, 

C 

K, 


atm 


flow area 

defined as intercept of linear velocity distribution in z 

defined as slope of linear velocity distribution in z 

2 

defined as 18*pg/d *p 
specific heat at constant pressure 
specific heat at constant volume 
nozzle diameter 
particle diameter 

friction factor for tube flow (=0.02) 
isentropic flow function (in /g /AP ) 

o c o 2 

gravitational constant (32.2 lb -f t/lb^-sec ) 

specific heat ratio (C /C ) 

P ^ 

particle to gas velocity ratio (V /V ) 

P 8 

particle to gas velocity ratio in nozzle section 
particle to gas velocity ratio in straight section 
length of straight section 

fictitious length necessary to achieve M g = 1. 
mass flow rate 

mass of particles per mass of gas (lb m particles/lb m gas) 
Mach number 
static pressure 
atmospheric pressure (=P^) 


volume flow rate 
gas constant 
temperature 
velocity 

axial distance from inlet 
Starred quantities - gas dynamics reference values 
e density 

u absolute viscosity 

u as unit microns 
Subscripts 
g - gas 
p - particles 
o - stagnation state 
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Subscripts continued 

i - inlet or stagnation conditions 
e - exit 

1 - upstream station 

2 - throat station 

3 - straight section exit station 
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I . INTRODUCTION 


The analyses and computations presented here are concerned with the 
description of the flow characteristics of two vacuum type particle 
accelerators which have been used in a series of experiments utilizing 
an MOS type impactor for atmospheric particulate detection. The accelerators 
are essentially vacuum driven nozzles which accelerate sample air and 
particulates to a high velocity. The high velocity mixture is subsequently 
stagnated on the sensitive portion of the MOS impactor, at which point each 
of the particles generates a discharge which is ultimately counted in the 
MOS circuitry. In order to understand the important mechanisms occurring 
in the impaction process and to properly interpret the results of the 
experiments, it was necessary to know the actual velocities of the particles 
at impact under experimental conditions. The task of the work presented 
here, therefore, was to analyze the two types of nozzles which were used 
under vacuum conditions and to make estimates of the actual gas and particle 
velocities which were attained under experimental conditions. 

Before gas dynamic type calculations could be carried out on the nozzles 
used in the experiments, it was necessary to investigate the accuracy which 
could be expected from a purely gas dynamic computation involving a particle 
laiden gas. This was done and is the subject of section II, p. 60. 

The first nozzle configuration used in the experiments was a converging 
nozzle, designed essentially to accelerate the gas-particle sample to a 
velocity sufficiently great to discharge the detector. A description of 
this device and some calculations of its characteristics is included in 
section III, p. 62. 

It was realized that this rapid acceleration of the gas-particle mixture 
tended to accelerate the smaller particles more than the larger ones, thus 
tending to eliminate the larger particles role in detection. A combination 
of converging nozzle and straight tube was then used to provide a slower 
acceleration over a longer distance. This configuration tends to accelerate 
the larger particles significantly more. An analysis of the gas dynamic 
characteristics of this configuration is included in Section IV, p. 66, and 
a particle dynamics analysis is included in Section V, p. 70. 


61 



II. EFFECTS OF PARTICLES ON GAS DYNAMICS 


The objective of the material presented here is to show that when 
particle concentrations are as low as they are in atmospheric air, gas 
velocities can be computed from gas dynamic computations by neglecting 
the presence of particles. 

For internal flows, the conservation of mass of the gas phase is 


m = p V A 
g g g 


( 1 ) 


and the conservation of linear momentum for the composite medium is 


m 


dV dv ~ 

l -r-®- + m. + A— = 0 

g dz p dz dz 


( 2 ) 


Assume that the ratio of the velocities of the particles and gas is 
K, i.e. 


K = V /V = const. 
P 8 


(3) 


This can be proven to be approximately true for nozzle flows (see 
ref. 1). Substituting equations (1) and (3) into equation (2) yields 


dV 

Og s g jf- a + V ) + f - 0 


M> 


where is (p /p ) and p^ is the "partial" density of the particles. 

Thus, m is lb particles/lb gas. Now, this is the one-dimensional 
p m r m 

momentum equation of gas dynamics when 


T 


g g 


p (1 + m K) 


(5) 


Also, p =pRT=p'R T/(Hm K) = p'R'T where 
g g g g g P g g 


R' = R / (1 + m K) 
g g P 


( 6 ) 
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Thus, it can be shown that if K = constant, the gas dynamic energy equation 
can be used in its homogeneous form when we replace unprimed by the following 
primed quantities: 


C' 

P 


C + m C (T -T )/(T -T ) 
= P§ P BE 2 P 2 § 

1 + m K 2 
P 


(7) 


= C' - R' 

(8) 

P 

= C 'J C l 

P D 

(9) 


These primed quantities are known as "pseudo" properties in two phase 
flow theory. Equations (5) through (9) show that when m^ is small the 
flow characteristics of a nozzle are the same as if the particles were 
absent altogether, an assumption which is valid for the relatively dilute 
particulate concentrations in atmospheric air. Thus, particle sampler 
flows can be treated first from a gas dynamic point of view and secondly 
from a particle dynamic point of view. This is the computation approach 
which has been used here. 
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III. VACUUM NOZZLE CHARACTERISTICS 


The purpose of the computations presented in this section was to 
first define the theoretical flow characteristics of the given nozzle 
and secondly to estimate actual gas velocities for the given nozzle 
with measured sampling rates. The computations make use of classical 
gas dynamic theory and procedures as presented in references 2,3,4. 

A description of the nozzle follows, along with the essentials of the 
compuations . 

The vacuum nozzle is a converging, conical nozzle with an inlet 

diameter of 0.225 inches, an exit diameter of 0.060 inches, and a 60° 

taper. A vacuum pump is used to establish a vacuum pressure at the 

exit area while the larger area is exposed to the sample gas at atmospheric 

pressure. The inlet area is assumed to be much larger than the exit area, 

2 

the latter of which is 0.00283 in . The stagnation density at the entrance 

3 

area is taken as 0.075 lb^/ft , the stagnation pressure at 14.7 psi and 
temperature of 70°F. The exit pressure, established by vacuum pump, and 
measured flow rates are as follows: 

P = 13.5 in. Hg. Vac.; 17.5 liter/min. 

P = 10 in. Hg. Vac.; 16.2 liter/min. 

P =2 in. Hg. Vac.; 7.8 liter/min. 
e 

A. Theoretical Flow Characteristics of Given Nozzle - In order to 
calculate the exit Mach number, M , the flow rate, Q, and the exit 
velocity, V , it was necessary to make the assumption that the nozzle 
behaves as an ideal (isentropic) expansion device from atmospheric 
conditions to the prescribed exit conditions. The pressure ratio P^/P^ 
is first computed and M^ determined from gas dynamic tables. The 
quantity m/RT^/ g^/AP^ is determined from Fig. 6-4 of reference 3, noting 

that P q = P^. The quantity /RT Q /g */A e P o was computed to be 712.0 for the 
experimental conditions (T -T^) , and m subsequently calculated. The 
density ratio p/P Q was determined from the gas dynamic tables and p, Q and V 
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then calculated. Note that the Q evaluated at both inlet and exit 
conditions is given. A summary of the computations are shown in the 
following table. 


P^Cin. Hg. Vac. ) 

P /P. 
e i 

M 

e 

m/RT /g '/AP. 

O Li 1 

m 

P/Po 

£ 

% 


V 

e 

13.5 

0.549 

0.97 

0.68 

9.55(10)~ 4 

0.65 

.0487 

33.3 

21.5 

998 

10 

0. 666 

0.785 

0.65 

9.13Q0)" 4 

0.75 

.0562 

27.5 

20.6 

823 

2 

0.933 

0.315 

0.34 

4. 77 (10)~ 4 

0.95 

.0712 

11.4 

10.8 

341 

The flow rates, 

Q are 

given 

in liter/min. 

the mass flow 

rate 

m is lb /sec 
m 

and 



the exit velocity given in ft/sec. The flow rate should correspond to 

the experimentally determined sampling rate. A plot of and versus 

P /P^ is shown in Figure A-l. Comparing the Q i (liters/min. ) with the measured 

values, it is seen that the theoretically predicted flow rate is approximately 

25% larger than measured. This was expected since inefficiencies in a simple 

nozzle such as that used in the experiments are normally of this order. 

B. Actual Velocities of Given Nozzle for Measured Flow Rates - The 

actual velocities at the nozzle exit section which correspond to actual 

measured inlet flow rates were computed. From the measured Q and the inlet 

density, m is computed. From this m and again the inlet conditions, the 

flow function m/RT /g '/A P is computed. This function is tabulated in 

o °C e o r 

reference 4 versus M . Thus, with this M the Q can be determined from 

e e e 

p/p Q and thus V^. Note that this procedure does not require an isentropic 
computation and therefore it is actual, then should be a good estimate. 

A tabulation of some computations are as follows. 



m 

t 

m/RT /g / AP 

O L* O 

M 

e 

e/e 

0 

^e 

V 

e 

17.5 

7. 75(10) 

0.552 

0.56 

0.86 

0.012 

610 

16.2 

7.18(10)" 4 

0.511 

0.50 

0.885 

0.0108 

550 

7.8 

3. 45(10)~ 4 

0.246 

0.21 

0.978 

0.0047 

239 
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3 

The flow rate Q. is in liters/min, m in lb /sec, Q in ft /sec and V 
1 me e 

in ft/sec. Note that the actual exit velocity is only about 60% of that 

predicted theoretically, again an expected result. A plot of the 

actual Q. and V are also included in Figure A-l. 
i e 
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(V&Q Theory) 


O L = 5.625" I 

I V&Q Experi- 

X L = 10.00" ( mental 



L = 5.625" 
f = 0.02 

Isentropic nozzle 


L= 10.00" 
f = 0.02 

Isentropic nozzle 



Figure A-l. Nozzle Flow and Velocity versus Nozzle Pressure 



IV. CHARACTERISTICS OF NOZZLE AND STRAIGHT TUBE COMBINATION 


It has been recognized that even though a simple converging nozzle 
is sufficient to accelerate the gas, it may not provide a long enough 
flow path to allow time for the acceleration of entrained particles. 

Thus, it was thought that an additional length of straight tube attached 
to the exit of a converging nozzle would facilitate the acceleration of a 
greater percentage of the larger particles. A sketch of this configuration 
is shown inFigureA-2. The nozzle exit. Station 2, and the straight tube 
were of an 0.062 in inside diameter. Two lengths were used in the 
experiments, 5.625 inch and a 10.0 inch section. Again, the nozzle inlet 
diameter was large in comparison with the throat diameter. As before, 
a theoretical analysis of this nozzle configuration and computations of 
actual characteristics were carried out. A summary of some experimental 
data is as follows : 



l 2 = 10.0 

in. 


L 2 = 5.625 

in. 

p 

e 

V?i 


P 

e 

P /P. 
e 1 


18.5 

0.382 

12.6 

19.5 

0.35 

13.25 

14 

0.533 

11.8 

15 

0.50 

12.6 

4 

0.867 

6.8 

10 

0. 666 

11.1 




4 

0.867 

7.2 


P g is in inches of Hg. , and in liters/min. sampled. 

A. Theoretical Flow Characteristics of Given Combination - Assuming 

isentropic flow in the nozzle section and friction flow in the straight 

section, the computation procedure is of a trial-and-error nature. The 

Mach number at the exit station, 3, M^ is first assumed and the friction 

length fL^/D read from the Fanno tables of reference 4. The actual friction 

length fL 2 /D is then computed and added to fL^/D to form the friction length 

corresponding to station 2, f (L^+L^) /D. Then M 2 is read from the Fanno tables 

for this friction length, along with P„/p*. Since is known, P 9 /P„ can 

P 1 J P P J 

be calculated from C^) / (P^/P*) and ^ 2 ^ a.t\a ca i cu i at ed from (~) (~ ). 

3 atm 
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From P 2 /P atm > is determined from the isentropic tables and compared with 

the previously computed M 2 . Using this procedure, the following data was 

generated. Noting also that m = AP f //kT / g = f /712 and the flow function 

f is read from the isentropic tables of reference 4 for each M,,. 
w K 2 


Case 

P,/P „ 

3 atm 

h 

fL 2 /D 

h 

h 

P 2 /P 3 

P„/P 
2 atm 

f 

w 

m(10) 3 

^i 

1 

0.35 

5.625 

1782 

1.0 

0.43 

2.50 

0.875 

0.464 

0.652 

14.8 

2 

0.50 

5.625 

1.82 

0.72 

0.42 

1.77 

0.885 

0.448 

0.63 

14.2 

3 

0.666 

5.625, 

1.82 

0.51 

0.375 

1.38 

0.913 

0.404 

0.566 

12.8 

4 

0.867 

5.625 

1.82 

0.23 

0.254 

1.10 

0.956 

0.29 

0.407 

9.18 

5 

0.382 

10.00 

3.24 

0.81 

0.354 

2.40 

0.916 

0.39 

0.548 

12.4 

6 

0.867 

10.00 

3.24 

0.57 

0.336 

1.73 

0.920 

0.376 

0.528 

12.0 

7 

0.867 

10.00 

3.24 

0.23 

0.205 

1.12 

0.975 

0.231 

0.325 

7.35 


In the above tabulation, m is in lb /sec and Q. = m/0.075 is converted 

m 1 

to liters/min. Some additional computations of velocities were made and are 
tabulated below: 


Case 

^1 l 

2^ p atm 


12 

v 2 /v* 

M 

_3 

v 3 /v* 

V /V 

y 2 

13 

1 

. 87 (10) -2 

.9094 

. 958(10) _2 

487 

.4729 

1.0 

1.0 

2.11 

1030 

2 

. 84 (10~ 2 

.9170 

. 916 (10) -2 

467 

.4522 

.72 

.7508 

1.66 

775 

3 

. 755 (10) -2 
-2 

.9347 

. 808(10)” 2 
-2 

411 

.3999 

.51 

.5447 

1.36 

560 

4 

.542(10) 

.968 

.56(10) 

285 

.2770 

.28 

.3043 

1.10 

314 

5 

. 73(10)~ 2 

.9395 

. 777 (10) -2 

396 

.3830 

.81 

.8343 

2.18 

863 

6 

. 705 (10) -2 
_2 

.9445 

. 747 (10) -2 
-2 

380 

.3682 

.57 

.6051 

1.64 

623 

7 

.433(10) 

.9803 

.442(10) 

225 

.2182 

.23 

.2506 

1.15 

259 

In this 

table Q 1 and Q are flow rates evaluated 

at stations 1 

and 2 


respectively in ft 

/ sec. 

The density 

ratio 

1 was taken from the 

isentropic 


tables at M 2 and V 2 computed from The velocity ratios with the 

starred reference velocities were taken from the Fanno tables at M 2 and 
M^ and subsequently calculatedin ft/sec. In all of the computations, 
a friction factor f of 0.02 was used and checked against Reynolds number 
as being representative. The results of the computations are illustrated 
in Figure A-l. 
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B. Actual Velocities of Given Combination with Measured Flow Rates 


The computations presented here are the results of the calculations 
of the actual exit velocities for certain measured flow rates. As above, a 
trial- and-error procedure is necessary. The governing equations are 


f 

w 


(P_ /PJ 

atm J 
/k 


= M 3 (l + 


k-1 

2 




( 10 ) 


V 3 - M 3 (l + if M 3 2 r 1/2 Ag c RT 1 


( 11 ) 


and 

f = m(712) (12) 

w 


From a measured Q and a ! 

known p , m is determined and 

subsequently 

f . Then, 
w 

f (P Qf . /v./A 

w atm 3 

is computed from 

measured pressures and M 

determined 

from equation 

(10) by trial-and- 

-error. Subsequently, 

is computed 

from equation 

(11) . A summary 

of the computations are 

included in 

the following table. 

, (P atn/ P 3 ) 



Case 

m 

w A 


h 

1 

. 587 (10) -3 

1.0 

.92 

960 

2 

. 558(10) 3 

0.672 

.65 

705 . 

3 

. 491 (10) -3 

0.443 

.44 

487 

4 

. 369 (10) -3 

0.256 

.25 

280 

5 

. 557 (10)~ 3 

0.879 

.82 

869 

6 

. 552 (10) -3 

0.623 

.60 

654 

7 

. 301(10) _3 

0.209 

.21 

235 


In this table m is in lb /sec. and V„ in ft/sec. These values are shown 

m 3 

plotted in Figure A-2 along with the theoretical results. The experimental 
is lower than the theoretical for = 5.625 in. but in essential 
agreement for = 10.00 in. This is also true for the plot. This is 
believed due to the reduction of losses in the nozzle reaction which occur 
when the velocities in the nozzle section are reduced by adding additional 
length on the straight section. Computed values of were not included 
as they can be adequately estimated from theoretical values shown. 
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V. PARTICLE VELOCITIES IN NOZZLE AND STRAIGHT TUBE COMBINATION 

It can be shown (see reference 1, p. 211) that the ratio of particle 
to gas velocity in a duct is a constant when the axial variation of the 
gas velocity is linear. That is, 


V f 


= K = constant 


(13) 


when V = a + bz where z is axial distance from duct inlet. This result, 
g 

of course, is valid for low particle Reynolds numbers and Stokesean drag. 
The velocity distribution in z in nozzle and straight duct flows is here 
approximated to be linear. Although this isn't entirely correct, it is 
a good approximation for most nozzles and ducts and especially for the 
problem under consideration here. Thus, 


V v = KV 
P g 


(14) 


and 


= C, 


1-K 

K 


dz + V p ( Zl ) 


v p ■ c i < ¥ i (2 - z i> + w 


(15) 


Equation (15) is developed from the particle momentum equation incorporating 

2 

Stokes drag law. C^ = 18*)jg/d ’p p where d is the particle diameter. 


A. Entrance Nozzle (K = IQ 


(1-K ) 

v p - c i~ icT - 2 + V 2 - 0) 


(16) 


where Vp(z=0) = 0 and V p = . Now, » using this latter ratio 

(1-K ) V (0) (1-K ) 

v s " c ! ~pr ~ z + ~q ~ " c i ~ 2 2 + Y 0) (17> 
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Now, assuming V = V (0) + bz, it can be seen that 
g S 


b = 


C 1 (1-K 2 ) 

K _ 2 


or, taking the positive root, 


K 2 


-C 1 + /C 1 + 4bC 1 

2b 


( 18 ) 


and 


V -V (0) 

b --8-B 

z-0 


(19) 


Knowing V g2 and z 2 , b can be determined and thus K 2> From , then 

V can be calculated. From the entrance nozzle V (0) = 0 and V (0) = 0. 

* C (1-K ) 6 “ 

Thus, b = V g2 /z 2 and V p = K 2 - 1 -~ 2 2 - z. 


Thus , 


K„ 


(i-k 9 ) 

v p 3 c i ~kT 


( 20 ) 


where C ;L = 18 yg/d p p and z 2 = 3/32 in. = 0.0078 ft. 


Computations were made for three particle sizes. These, along with some 
property data are listed below. 


d = 3p = 1.18(10)" 4 in. 

=9 y = 3.55(10)' 4 in. 

= 15y = 5. 91 (10) _4 in. 

p p = 2.5(62.4) = 210 lb m /ft 3 

V . = 0.042 lb /hr. ft. 

air m _ 

C = 1.38(10)“ 4 /d sec 


d 

1. 18(10)~ 4 

2.36(10) 4 

3.55(10) 4 

5.9K10)' 4 

in. 

C 1 

0.993(10) 4 

. 248(10) 4 

.1095 (10) 4 

0.0395(10) 4 

-1 

sec 


The computations are summarized in section C, p. 73. 
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B. Straight Tube (K=K J 


A similar type of analysis is carried out for the straight tube. 
However, for the tube is larger in magnitude than K 2 for the nozzle. 
Thus, at station 2 V can have two values depending on the value of K 
used to compute it. A correction is applied to avoid this difficulty. 
The particle velocity is given by 

<1-K,) 

Vp - c rir— < z - z 2> + W < 21 > 


and the gas velocity by 


(l-K,) 

V g - Cj— j— < 2 -b 2 ) + V p(22 )/K 3 
K „ 


( 22 ) 


Again V = V (z.) + b(z-z„) where b = (V _-V „)/z„-z»). Thus, 
g g 2 2 g-> g 2 3 2 

(1-K ) 

V g (z 2 ) + b(z-z 2 ) = C ± y~ (z-z 2 ) + V p (z 2 )/K 3 . 


(23) 


Now 


, if V (z 2 ) is placed equal to V p (z 2 )/K 3 then as before 


bK„ = C 


(1-K 3 ) 


3 1 K. 


This however would produce a V p (z 2 ) which is too large because is 
larger than K 2> Also, V p (z 3 ) would be too large by the same amount 
because the distribution is linear. Fortunately, this excess in V p in 
the straight section is known and can be subtracted from the final result. 
Thus 


V = K *V - (K -K )V 
P3 3 g3 k 3 2 J g2 


(24) 


where b^ 


(Vg3 V g2 )/(z 3 z 2^ anc * 


K 3 


- C 1 + /C 1 + 4b 3 C 1 

— 


(25) 


The straight tube computations are also summarized in the following section. 
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C. Computations of Particle Velocities 

The following is a summary of the computations outlined in 
sections A and B. 


d~3p 


V p i 

12 

fa 

(z 2 -°) 

(z 3 -z 2 ) 

fa 

fa 

fl fa 

fa 

V P 2 

V P 3 

0.35 

480 

960 

0.094 

5.625 

1025 

6 . 12 ( 10) 4 

• 993 ( 10) 4 .912 

.33 

157 

596 

0.50 

450 

705 

0.094 

5.625 

545 

5 . 75 ( 10) 4 

. 993 ( 10) 4 .945 

.338 

152 

392 

0.666 

375 

487 

0.094 

5.625 

239 

4. 80 ( 10) 4 

. 993 ( 10) 4 .985 

.364 

136 

245 

0.867 

260 

280 

0.094 

5.625 

42.5 

3 . 32 ( 10) 4 

. 993 ( 10) 4 1.0 

.417 

109 

128 

0.382 

400 

869 

0.094 

10.00 

563 

5 . 1 ( 10) 4 

. 993 ( 10) 4 .942 

.354 

141 

585 

0.533 

380 

654 

0.094 

10.00 

329 

4 . 85 ( 10) 4 

. 993 ( 10) 4 .973 

.362 

137 

404 

0.867 

200 

235 

0.094 

10.00 

42 

2 . 55 ( 10) 4 

• 

. 993 ( 10) 4 1.0 

.47 

88 

129 


d = 9 p 





d = 15 y 




P /P 

y 1 

fl 

K 3 K 2 

V P2 

V P 3 

fi fa 

f 2 

V P 2 

V P 3 

.35 

. 1095 ( 10) 4 

.63 .125 

60 

363 

3 . 95 ( 10 ) 2 .454 

.0771 

37 

255 

.50 

. 1095 ( 10) 4 

.732 .129 

57.8 

244 

3 . 95 ( 10) 2 .56 

.0795 

35.7 

179 

.666 

. 1095 ( 10) 4 

.837 .14 

52.5 

147 

3 . 95 ( 10) 2 .693 

.0867 

32.5 

110 

.867 

. 1095 ( 10) 4 

.954 .167 

43.5 

62 

3 . 95 ( 10) 2 .856 

.103 

26.8 

0 

.382 

. 1095 ( 10) 4 

.728 .146 

58.4 

399 

3 . 95 ( 10) 2 .554 

.0842 

33.7 

294 

.533 

. 1095 ( 10) 4 

.807 .142 

53.9 

274 

3 . 95 ( 10) 2 .641 

.0862 

32.7 

208 

.867 

. 1095 ( 10) 4 

.893 .188 

37.6 

79 

3. 95 ( 10) 2 .845 

.117 

23.4 

52 


In the above tabulation V and V„ are the gas velocities in ft/sec., the 

^ ^ —1 —1 

z's are in inches, b^ and b 2 in sec , in sec and the particle velocities 
in ft/sec. These latter quantities were computed by means of equations (24) 
and (20) and are plotted in Figure A-3. 

The accuracy of the above calculations is difficult to assess. Predicting 
particle entrainment in gas streams is a difficult problem. It is believed 
that the above computations at least establish a lower bound for the exit 
particle velocities. In all probability, the actual velocities will be of 
the order of the above estimates or slightly greater. The scheme used in 
the above computations should at least be accurate than a prediction based 
on a constant, average gas velocity, a computation which might be considered 
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Ft/sec 







a first estimate procedure. A more sophisticated procedure than that used 
above for making the computations is certainly not evident. 
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